Packing-limited growth of irregular objects 
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We study growth limited by packing for irregular objects in two dimensions. We generate packings 
by seeding objects randomly in time and space and allowing each object to grow until it collides with 
another object. The objects we consider allow us to investigate the separate effects of anisotropy 
and non-unit aspect ratio. By means of a connection to the decay of pore-space volume, we measure 
power law exponents for the object size distribution. We carry out a scaling analysis, showing that 
it provides an upper bound for the size distribution exponent. We find that while the details of the 
growth mechanism are irrelevant, the exponent is strongly shape dependent. Potential applications 
lie in ecological and biological environments where sessile organisms compete for limited space as 
they grow. 

PACS numbers: 02.70.Rr, 05.10.-a, 87.23.Cc, 87.23.-n, 81.10.Aj 
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I. INTRODUCTION 

In a previous work J5J], we examined packings formed 
by spheres growing in d dimensions, immediately stop- 
ping upon contact with another sphere. We termed this 
to be packing-limited growth (PLG). Here, we address 
the question of what happens when we consider irregular 
(i.e., non-spherical) objects of similar shape growing in 
d = 2 dimensions. We focus in particular on collisions of 
growing rectangles. By combining rectangles, we are able 
to form and examine a range of shapes whose growth pat- 
terns vary broadly in radial anisotropy. Whereas in [l]] 
we found universality classes depending only on dimen- 
sion d, we find here that marked non-universal behavior 
arises when shapes are varied. 

Our problem finds several motivations. First, this is a 
relatively unexplored kind of packing. Packings are gen- 
erally static or randomly generated [j2|: physical mecha- 
nisms are only seldomly connected with the creation of 
packings || Q. Furthermore, packings are typically of 
monodisperse objects of the same form or taken from a 
small set of forms [| . From a physical point of view, 
two-dimensional packings of growing objects may be of 
use in modeling or understanding certain biological and 
ecological patterns. We consider the geometric approach 
presented here as a preliminary step towards describing 
how shape alters the size distributions of populations. 

In section |ll| we describe the various objects we con- 
struct from rectangles. Although not all of these objects 
have obvious physical parallels, they present a range of 
typological cases from which direct app lications and com- 
parisons may be sought. In section III, we provide a scal- 
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ing analysis that empirically appears to be exact for d > 4 
in the case of hyperspheres H (details of the calculations 
for this section are given in Appendix [X]). We report 
the results of our numerical investigations in section IV 
along with a discussion of the failings of the scaling the- 
ory. We conclude the paper in section ^ and outline our 
algorithm for packing rectangles in Appendix M. 



II. IRREGULAR OBJECTS CONSIDERED 

The building block shape of the objects we consider is 
the rectangle. We denote the aspect ratio as a, defining 
all rectangles such that a > 1. We record rectangle size 
r as half the length of the long side, the dimensions then 
being 2r and 2r/a. Examples of each object are provided 
in the packings of Figures [I] and ||. 

For a = 1, we have squares and as a — > oo, rectan- 
gles effectively become line segments and the packing 
becomes one of fitting d = 1 objects into a d — 2 vol- 
ume. Rectangles afford a basic example of anisotropic 
growth since the long side grows at a rate slower than 
the short side by a factor of 1/a. The growing square is 
distinguished from the rectangle since the former expands 
uniformly perpendicular to its edges. However, in com- 
parison to disks, the growth rate of edge points relative to 
the center of both squares and rectangles is non-uniform. 

The simplest combination of two rectangles is a cross 
(Figures |2|A and ||B) and a natural generalization is the 
2n-spoke object (with n — 1 being a rectangle and n = 2 
being a cross, see Figure |]C). Each 2n-spoke object is 
then part of a family of shapes indexed by a. 

The last object we consider is an eight-pointed star as 
shown in Figure |^D. This object is formed by two squares 
overlaid at angle of 7r/4 to each other. 

All objects are packed using the approach of Manna [Q : 
objects are added sequentially and allowed to instantly 
grow so as to just reach the existing packing structure. 
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FIG. 1: Rectangle packings of the unit square for varying 
aspect ratio a created using the packing-limited growth mech- 
anism. Each packing consists of 1000 rectangles with periodic 
boundary conditions being imposed. The aspect ratios cor- 
responding to the packings are (A) a = 1, (B) a — 3, (C) 
a = 10, and (D) a = 100. The first three packings are initial- 
ized with four randomly placed and oriented rectangles with 
longest side length 0.25, with eight such rectangles being used 
in the fourth packing. 



FIG. 2: Packings of irregular objects formed by combinations 
of rectangles. As per Figure [j], each packing contains 1000 
objects and the boundary conditions are periodic. In (A) 
and (B), the crosses are formed by two rectangles at right 
angles with aspect ratios a — 3 and a = 10 respectively. Plot 
(C) shows a packing of six-legged stars composed of three 
rectangles with aspect ratio a = 10. In (D), the packing object 
is an eight-pointed star formed by two overlapping squares set 
at an angle of n/4 to each other. 



We have observed and argued ffl that the value of the size 
distribution exponent and all other related exponents is 
independent of the growth dynamic. 

In the context of packing, spheres make for straight- 
forward calculations since the contact point between any 
two colliding spheres occurs along the line through their 
centers. Rectangles constitute a relatively simple gener- 
alization of spheres from a numerical point of view, hence 
our use of them here. 



III. SCALING THEORY 

In , we argued that a simple scaling assumption may 
be made regarding the form of P(r;n), the distribution 
of sphere sizes after n objects have been packed. We take 
P(r;n) to be described by a power law for radii above a 
cutoff value r r and uniform below: 



P(r; n) 



q-i r -i 



for r < r r 
—rJ 1 (-) a for r > r c 



(i) 



The tail of the distribution is fixed and the distribution 
fills in (r c decreases with n) with the (n + l)th object 
being chosen from within the pore space. The size distri- 
bution P(r; n) is connected to the probability of inserting 



an object of size r after n objects have been deposited, 
Pns{r,n). Since the probability of adding a sphere of 
vanishing radii must be proportional to the total surface 
area of the existing spheres S(n), we were able to esti- 

as m 



mate Pi ns (r; n 



Pins(r;n) = 



$(n) 



< r < r c , 



(2) 



where $(n) is the pore space volume. Using equation (|l|) 
to calculate S(n) and $(n) and requiring that Pi n8 (r;n) 
be normalized, we obtained the estimate a — 11/4 for 
disks. 

More generally in d = 2 dimensions, if we now write 
area as A — pr 2 and perimeter as L — or, we find the 
scaling exponent to be 



9cr + 4p 
3cr + 2p' 



(3) 



Setting (7 — 27r and p = n, we recover a = 11/4 for 
disks. Note that the probability of adding a non-spherical 
object is still proportional to S(n) in the limit of r — > 
0, regardless of the object's shape. For example, it is 
irrelevant that an added square of side length 2r may be 
oriented such that its center is between r and y/2r away 
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FIG. 3: For rectangle packings of 10 5 objects, plots of 3>(n), 
the decay of pore space as a function of number of rectangles 
added n, and N(r), the frequency of objects of size r, for 
rectangles. In both plots, the rectangles have aspect ratios 
a = 1 (solid line), a = 10 (dashed line) and a = 100 (dot-dash 
line). N(r) is binned in log-space for clarity. The exponent 
a (in N(r) oc r~ a ) increases with aspect ratio a, limiting to 
3 as a —* oo. Correspondingly, /3 decreases towards (see 
equation (Eh). 



from its point of contact with the existing packing. Also 
note that as a — + oo, the theoretical prediction tends to 
a = 3, in agreement with a previous determination of an 
upper bound for polydisperse packings || . Details of the 
calculations finding a for the objects we consider here are 
to be found in Appendix [X| 

All other exponents depend on a via simple scaling 
relations Q. In particular, the pore volume $(n) decays 
as n _/3 with the connection between (3 and a being 



We use this equation to calculate theoretical estimates of 
[3, denoting them by (3. Since direct measurement of (3 
is a significantly more robust exercise than determining 
a from P(r), we employ equation (EJ) in estimating a in 
the following section. 

IV. NUMERICAL RESULTS 

For each shape, we generate statistics for single pack- 
ings with 10 5 objects. Some example distributions taken 



object 






(3 


$ 


a 


& 


disk 






0.278(1) 


0.1429 


2.564(1) 


2.750 


square 






0.223(2) 


0.1429 


2.635(2) 


2.750 


rectangle 


a — 




0.207(2) 


0.1000 


2.656(2) 


2.818 


rectangle 


a = 


51 


0.145(2) 


0.0526 


2.746(2) 


2.900 


rectangle 


a — 


10) 


0.094(1) 


0.0294 


2.828(2) 


2.943 


rectangle 


a — 


20) 


0.055(1) 


0.0156 


2.897(2) 


2.969 


rectangle 


[a = 


•50) 


0.0242(3) 


0.0065 


2.953(1) 


2.987 


rectangle 


[a = 


100) 


0.0125(1) 


0.0033 


2.975(1) 


2.993 


cross (a = 


3) 




169 (3) 


0.0847 


2.710(3) 


2.844 


cross (a = 


10) 




0.074(1) 


0.0307 


2.862(2) 


2.940 


cross (a = 


100) 


0.009(1) 


0.0033 


2.983(2) 


2.993 


six spoke (a = 


10) 


0.078(2) 


0.0318 


2.855(3) 


2.938 


eight-pointed star 


0.213(2) 


0.1429 


2.648(4) 


2.750 



TABLE I: Estimates of a, the exponent of the number distri- 
bution, P(r) oc r~ a , for irregular objects undergoing packing- 
limited growth. The results for disks are included for compar- 
ison EJ. All other objects are combinations of rectangles (see 
Figures with a being the aspect ratio. The exponent 

/3 is determined from $(n) and a is subsequently obtained 
using equation (EJ). Full details of the method of measuring 
/3 are given in the text. Each measurement is for a single 
packing containing 10 objects. Measurement errors reflect 
variation in the approach to a limiting value of the exponents 
using the method of Section The scaling theory estimates 
of p and a, which are lower and upper bounds respective- 
ly, are calculated using equations (^|), (^), and those given in 
Appendix lAl 



from rectangle packings for a = 1, 10, and 100 are shown 
in Figure |3|. The distribution N(r) is the unnormalized 
frequency distribution corresponding to P(r). In the case 
of plain rectangles, we see that as a increases, both the 
observed a and the theoretical estimate a increase while 
(3 and (3 accordingly decrease. Furthermore, for all the 
shapes considered we find that a < a for all finite a. 

We note that measuring power law exponents is not a 
trivial procedure and here we take some care to ensure 
the validity of our results. While simple regression is 
the basic tool of analysis, the presence of, for exam- 
ple, crossovers and finite size cutoffs can substantially 
degrade the level of precision attainable. The method of 
measuring exponents we use here is based on examining 
a smoothed version of the derivative of the distribution 
as viewed in double logarithmic space (similar approach- 
es are to be found in [pj, 0). For example, for <&(n), we 
perform regression analysis on log 10 $(n) versus log 10 n 
over a sliding, variable width window of values of log 10 n. 
Writing the upper and lower limits of this window as 
log 10 ni and log 10 ri2, we have w = log 10 n^/rii being the 
width. In general, depending on the number of orders 
of magnitude spanned by the data, we would preferably 
choose w in the range 0.5 < w < 3. Here, we fix w — 1 
and calculate the "local" exponent (3(ni) for each win- 
dow. We find in all cases that (3{n\) tends towards a 
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log 1Q a 



to be independent of growth dynamics. However, the 
main result presented here is that the exponent is highly 
shape-dependent, adopting a continuous range of values 
2.564 < a < 3. Ultimately, understanding how geom- 
etry impacts on the structure of plant communities will 
require imposing a notion of packing-limited growth onto 
a reasonable set of dynamics. The results here demon- 
strate that in so doing we must keep in mind that shape 
matters. 
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FIG. 4: Measured values of a as a function of aspect ratio a 
and shape. The symbols correspond to disks (circles), rectan- 
gles (squares), crosses (plus sign), six-spoke (hexagon), and 
eight-pointed star (asterisk). 



constant value. This indicates the scaling law is robust 
and further allows us to estimate (3 along with an error 
based on the fluctuations observed for j3(n\). Measured 
values of a and (3 are recorded in Table |. 

As we have noted above, our scaling theory approach 
appears to be exact for d > 4 dimensions in the case 
of hyperspheres and an overestimate of the true value 
of a for d < 4. We therefore do not expect the scaling 
theory to be correct for non-spherical objects in d = 2 
dimensions. Indeed, in all cases, we observe the theo- 
retical estimate a is an overestimate of the measured a. 
This direction of error makes sense in light of the scaling 
assumption made by equation (||). The actual form of 
Pins(r]n) is not precisely uniform but rather rolls over 
less steeply than a step function around r = r c . Given 
the manipulations that leads to (||) , it can be argued that 
a > a must hold Q. 

In Figure we show how a varies as a function of 
both aspect ratio a of the constituent rectangles and the 
object's particular shape. The strongest influence is evi- 
dently aspect ratio with a varying from 2.564 (disks, 
a = 1) up to 3 (rectangles, n-spokes, a — oo). We 
observe a secondary effect due to the details of the shape. 
Square and the eight-pointed star packings have a value 
of a increased above that of disks. Cross packings have a 
higher value of a than do rectangles with the same aspect 
ratio. 



V. CONCLUDING REMARKS 

We have extended a model describing the interaction 
of growing disks Jlj] to the problem of growing irregu- 
lar and anisotropic objects. In this model, the exponent 
characterizing the size distribution of objects is found 



APPENDIX A: SHAPE PARAMETERS FOR 
MEAN FIELD CALCULATIONS 

In this section, we derive the formulas for the area 
and perimeter of irregular objects in d = 2. For general 
rectangles, we have area given by 

A-cct = 4r 2 /a, (Al) 

and perimeter as 

ircct = 4(1 + l/a)r. (A2) 

Using equation (|^) , we therefore have the mean field esti- 
mate of a for rectangles as 

Ckroct = 3 — 



5 + 3a' 



(A3) 



where writing the result in this fashion makes plain the 
limiting value of d re ct = 3 for a — > oo . For the 2n-spoke 
objects described above, area grows as 

-4n-spoko = 2n/a 2 (2a - cot7r/2n)r 2 , 



and perimeter is given by 

in- S poko = 4rt/a(a + 1 - cot7r/2n)r, 
leading to the estimate 

_ 13 + 9a - (9 + 2/a) cot 7r/2n 

^n-spokc 



(A4) 
(A5) 

(A6) 



5 + 3a- (3 + l/a)cot7r/2n ' 

Note that fo r n = 1, we recover the result for rectangles 
in equation (A3). Finally, for the eight-pointed star we 



have 



and 



A 8 - P oint = 2[2 + (l-tan7r/8) 2 ] 



is-point = 16(1 - tan7r/8)r, 



which yields 



42- 40 tan ttJ 



2tan 2 7r/ 



^8-point 



15 - 14tan7r/8 + t air 717 
the same as for disks and squares. 



(A7) 
(A8) 

11/8, (A9) 
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APPENDIX B: RECTANGLE COLLISIONS 

We describe an arbitrary rectangle in the x-y plane as 
follows. The sides of the rectangle are 2r and 2r/a where 
r > is half the length of the long side of the rectangle 
and a > 1 is take the aspect ratio. The rectangle is 
centered at (xq, yo) an d rotated at an angle 8. We take 8 
as the angle between the direction of the positive x-axis 
and the short axis of the rectangle (i.e., parallel to the 
side with length 2r/a) so that < 8 < it. 

A rectangle described by (r, a) and 9 — and (xo, yo) — 
(0, 0) satisfies the equation 



max 



ax 
I) 



a 
1 



(Bl) 



This is a remapping of a rectangle into a square of side 
length 2r. For an arbitrary rectangle, we can map it onto 
a basic square by recentering it at the origin, removing 
the rotation, and undoing the dilation by the aspect ratio 
a: 



max 



a 
1 



cos 9 sin 9 
— sin 9 cos 9 



x — xo 

y-yo 



r. (B2) 



Therefore, to determine whether or not an arbitrary 
point (x, y) lies on or within a given rectangle, we need 
to check whether or not 



a 
1 



cos 9 sin f 
— sin 9 cos I 



X — Xq 

y-yo 



< r. (B3) 



In packing growing rectangles, the above is used to check 
if a newly seeded rectangle has been placed in pore space 
and not within an existing rectangle. 

For a rectangle that passes this test, the next calcu- 
lation is to determine how large it may grow preserving 
its aspect ratio and orientation so that it just reaches an 
existing rectangle. 

To do this, we consider one arbitrary rectangle that 
does not cover the origin and find the size of largest rect- 
angle centered at the origin with the short side along the 
x-axis such that the rectangles just touch. We will then 
generalize to any configuration by appropriate rotations. 

First, the "growing" rectangle could hit the already 
existing one at any of the latter's corner points. The 
four corner points are given by 



x 


+ r 


. y° . 





x 




+ r 


. y° . 





± 1/a cos 9 — sin 9 
± 1/a sin + cos 9 

=Fa cos 9 + sin 9 
=Fa sin 9 — cos 9 



(B4) 



(B5) 



The other possible collisions are between the corners of 
the added rectangle and the sides of the existing rectan- 
gle. There may be 0, 1, 2, 3 or 4 such interceptions (1 or 
3 if the rectangles touch at corners). To calculate these 



points, we parametrize each side of the existing rectan- 
gle and find the interceptions with y — ±ax. For an 
arbitrary line segment described by 



x = bi + b^t and y = c\ 



(B6) 



with — 1 < t < 1, the intersection with y = ±ax occurs 
when 



t 



— ci ±abi 



(B7) 



Parametrizing the sides of an arbitrary rectangle gives 
the following: 



X 




Xq 


+ r 


. y . 








X 




x 


+ r 


. y . 








X 




X 


r 

+ - 


y 




yo 


a 



— sin# 
cos 9 

sin (9 

- cos 9 



rt 
a 

rt 
a 



cos 9 


+ rt 




sin 9 





cos 9 
sin6> 

cos 9 
sin 8 

— sin 9 
cos 8 



, (B8) 
, (B9) 
, (BIO) 
.(Bll) 



The first two equations describe the short sides and the 
latter two describe the long ones. 

Using equation (B7) in equations (B8) to ( Bll| ), we 
respectively have eight possible solutions (2 for each 
side) : 



X 




x 


r 

+ - 


— cos 9 


+ rt 


— sin9 


. y . 




. yo . 


a 


— sin 9 




cos 9 



t = 
t = 



'(yo + rcos9) ± a(xo — rsin9) 

r/asin9 =F rcos9 
- (yo — r cos 9) ± a(xo + r sin 9) 
r/asm.9 =p rcos9 
-l/a(yo + r/asin9) ± (xo + r/acosff) 

r I a cos 9 =F r sin 9 
-\/a(yo — r/a sin 9) ± (xo — r/acos9) 
r/acos9 =F rsin(9 



(B12) 
(B13) 
(B14) 
(B15) 



Each of these has to be tested to see if —1 < t < 1. 
If so, then upon sub stituting t he va lues of t determ ined 
in in equ ation s (B12) through (B15) into equations (BS) 
through ( Bll ), we have 



(xo sin 9 — yo cos 9 — r) 

sin 9 =p a cos 9 
(xq sin 9 — yo cos 9 + r) 
sin 9 =F a cos 9 
(xq cos6> + yo sin9 + r/a) 

cos 9 ± a sin 8 
(xo cos 9 + yo sind — r/a) 
cos 8 ± a sin 8 



(B16) 
(B17) 
(B18) 
(B19) 



where y = ±ax in all cases. 

We are now able to write down how to determine the 
size of the largest rectangle that may fit in given it is 
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centered at the origin with an angle 4> and there is one 
other rectangle already present at (xq,uq) oriented at an 
angle 9. 

1. Determine whether or not the ne w r ectangle lies 
within the old one using equation (|B3|). 



2. If it is not enclosed, rotate the coordinate system 
by an angle <f> so the central rectangle sits square 
with the axes. All that needs to be done is to move 
the existing rectangle from (xo,yo) to (cos^xo + 
sin <fiy a , — sin (f>x + cos (f>y ), and change its angle 
of rotation to 6' = 9 — <\> (mod 7r). 



3. 



Find the pos ition of the four corner points using 
equation (B4) with 9' and (x' ,y' ). 



the diagonals of the new rectangle e xist by calcu- 
lating t in equations ( B12 ) through ( B15 ). 



5. Each value of t that satisfies — 1 <t < 1, is a valid 
intersection point. The positio ns of all such point s 
are then taken from equations ( |B16| ) through ( |B19| ) 
along with y = ±ax. 

6. For each valid (x,y) pair, calculate max(|ax|, \y\) 
which gives us the half-width of the shortest side 
of the rectangle centered at the origin that passes 
through the point (x,y). Take the minimum over 
all such values to find the half-width of the short- 
est side of the actual largest rectangle that may be 
inserted without overlapping. 



4. Determine which, if any, intersection points with 
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